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Abstract 



Various higher-dimensional black holes have been shown to be unstable by studying 

linearized gravitational perturbations. A simpler method for demonstrating instability 

is to find initial data that describes a small perturbation of the black hole and violates 

CJ ! a Penrose inequality. An easy way to construct initial data is by conformal rescaling of 

the unperturbed black hole initial data. For a compactified black string, we construct 

initial data which violates the inequality almost exactly where the Gregory-Laflamme 

instability appears. We then use the method to confirm the existence of the "ultraspin- 

^ ^ ning" instability of Myers-Perry black holes. Finally we study black rings. We show 

^ ■ that "fat" black rings are unstable. We find no evidence of any rotationally symmetric 

•^ , instability of "thin" black rings. 
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1 Introduction 



Vacuum black holes in four spacetime dimensions are believed to be stable against grav- 
itational perturbations. A qualitatively new feature that emerges in higher dimensions is 
the possibility of black objects with unstable horizons. The first example to be discovered 
r> I was the Gregory-Laflamme instability of a black string [1]. Later, heuristic arguments were 

c^ ' presented which suggest that " ultrasp inning" Myers- Perry [2] black holes should suffer from a 

similar kind of instability [3]. The existence of this instability has been confirmed by studies 
of linearized perturbations [U El El [7] . 

The equations governing linearized perturbations of higher-dimensional rotating black 
holes are very complicated. It would be nice if there were a simpler method of demonstrating 
black hole instabilities. In this paper, we will show that the existence of certain types of 
instability can be demonstrated using inequalities analogous to the Penrose inequality (see 
Ref. [8] for a review). In the form presented in Ref. |9], the 4d Penrose inequality is 

Amin < IGttE^ (1) 

where E is the ADM energy of asymptotically fiat initial data for Einstein's equation and Amin 
is the greatest lower bound on the area of any surface that encloses the apparent horizon of 
this data (Amin may be less than the area of the apparent horizon [lO]). This inequality results 
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from the assumption that the spacetime resulting from initial data containing an apparent 
horizon must "settle down" to the Kerr solution. It has been proved only for the special case 
of time-symmetric initial data [HI WI\ . 

The inequality ([1]) is supposed to apply to any asymptotically flat initial data containing 
an apparent horizon. One can also consider a less general local Penrose inequality, in which ([1]) 
is restricted to initial data describing a small perturbation of a Schwarzschild black hole. Such 
an inequality was first investigated by Gibbons [13], who considered certain time-symmetric 
initial data describing a small perturbation of the Schwarzschild solution. He found that ([T]) 
held to second order in perturbation theory for his data (at first order, the result follows from 
the first law of black hole mechanics [H]). Related results were obtained in Ref. [15]. 

In this paper, we will use local Penrose inequalities to demonstrate instability of certain 
stationary black hole solutions. If such a black hole is stable, then initial data describing a 
small perturbation of the black hole must satisfy a local Penrose inequality analogous to ([T]). 
Therefore, if we can find initial data that violates the inequality then the black hole must be 
unstable. 

Our first example of this method is the Gregory-Laflamme instability. Consider a Schwarzschild 
black string compactified on a Kaluza-Klein circle of radius 27tL. A surface of constant t in 
this geometry gives initial data for Einstein's equation. Now consider new initial data cor- 
responding to a small (but finite) perturbation of the black string (preserving the radius of 
the KK circle at infinity). Assume that the black string is stable. Then the perturbation 
will disperse through radiation to infinity and across the horizon so the evolution of this data 
will settle down to a new stationary black string. In general, this might have a small angular 
momentum Jp and a small linear momentum Pp. The solution will be a boosted Kerr black 
string with mass Mp and horizon area 

^F < -^ (2) 

The inequality follows from the fact that a boosted Kerr string has smaller horizon area than a 
Schwarzschild black string of the same mass. Let Aj denote the area of the intersection of the 
event horizon with the initial surface. Since the event horizon encloses the apparent horizon 
we have A^i^ < Aj < Ap where A-ann was defined above and the second inequality follows 
from the second law. Gravitational waves carry away energy (the Bondi energy decreases) so 
Mp < E where E is the ADM energy of the initial data. Combining these inequalities we are 
led to 

The inequality (E]) is a local Penrose inequality for the Schwarzschild black string. 

The assumption made in deriving ([3]) is that the string is stable. Hence if we can find 
suitable initial data which violates the inequality then we have demonstrated that the string 
cannot be stable. We shall construct initial data describing a small perturbation of the black 
string simply by conformal rescaling of the initial data for the unperturbed string. Our 
initial data violates the local Penrose inequality when r_|_/L is smaller than a certain critical 
value (where r_|_ is the horizon radius of the unperturbed string). The critical value for r_|_/L 
is smaller than the critical value at which the GL instability is known to appear by less 



than 0.2%. Hence the existence of the Gregory-Laflamme instabihty can be predicted from 
initial data alone. The close agreement between our result and the GL result is somewhat 
surprising: violation of the local Penrose inequality is a sufficient, but not necessary, condition 
for instability. 

Note that initial data corresponding to a constant t surface in the black string spacetime 
saturates the inequality (|3]). Hence, if the string is stable, such initial data minimizes E for 
fixed Amin or maximizes Amm for fixed E. In other words, a stable string must be a local 
minimum of energy for fixed horizon area or a local maximum of horizon area for fixed energy 
in the space of asymptotically KK initial data for Einstein's equation. Our proof of instability 
amounts to showing that certain black strings fail to satisfy these propertieso 

The main aim of this paper is to use the above argument to demonstrate instabilities of 
higher-dimensional rotating black holes. In the rotating case, one needs to assume that the 
initial data preserves some rotational symmetry in order to derive a useful Penrose inequality 
|13j . Therefore this method can be used to demonstrate instabilities that preserve some 
rotational symmetry. This includes the Myers- Perry instabilities discovered in Refs. [ll[5l[6],[7]. 
In this case, we shall confirm the existence of the "ultraspinning" instability predicted in Ref. 
[3] with far less effort than required for the numerical analyses of Refs. [H |H]. However 
we cannot predict instabilities which break the relevant rotational symmetry, such as the 
"bar- mode" Myers- Perry instabilities found in Refs. [TBI [T7] . 

Next we consider the stability of the black ring solution of Ref. [18] . So far, investigations 
of black ring stability have been heuristic. For a given mass, there is a finite range of angular 
momenta for which there exist two distinct ring solutions, referred to as "thin" and "fat" 
because of the shape of the horizon. Heuristic arguments suggest that "fat" rings should be 
unstable, as we now explain. 

Ref. [19] used "turning point" methods to argue that fat rings must have one more 
"unstable mode" than thin rings. This argument assumes that the "states or configurations 
of a given system" correspond to points in some manifold Ai on which one can define quantities 
such as the mass, angular momenta, and entropy. "Equilibrium states" correspond to points 
which extremize the entropy at fixed mass and angular momenta. These are assumed to form 
a submanifold A^cq of A^. In the present case, this corresponds to the known black ring 
solutions, hence one knows the entropy, etc., on A^eq- A given equilibrium state is stable if it 
is a local maximum of entropy for fixed mass and angular momentum. Stability can change at 
either a "turning point" or a "bifurcation". The former can be identified from knowledge of 
the thermodyamic quantities on A^cq- For black rings, there is a turning point as one moves 
from the thin ring to the fat ring branch. From this, one can deduce that fat rings near to 
the turning point are unstable. 

This argument is very suggestive (and successfully used in astrophysics to predict neutron 
star instabilities) but not entirely rigorous since, for example, the manifold M. is not defined. 
A related point is that it is not clear whether the predicted instability should be present 
classically. Finally, it only allows one to deduce that fat rings are unstable near the turning 
point. It cannot be concluded that all fat rings are unstable because it is possible that there 



^A similar argument was used in Ref. jl4| to explain why coloured black hole solutions of Einstein- Yang- 
Mills theory should be unstable. 



is another change in stabihty (at a bifurcation) as one moves along the fat ring branch. 

A different approach was taken in Ref. [2U], which considered certain singular deformations 
of the black ring to deduce an "effective potential" for variations of the radius of the ring. It 
was found that fat rings sit at a maximum of the potential and therefore should be unstable. 
Again, this is very suggestive (and in agreement with the result of Ref. [19]) but it is not clear 
that such a simple mechanical picture captures all of the relevant gravitational dynamics. 

We will use violation of a local Penrose inequality to show that all fat black rings are 
indeed classically unstable. We will also construct multi-parameter families of initial data 
describing perturbations of thin black rings. These all respect the local Penrose inequality 
and so our results are consistent with thin rings being stable against rotationally symmetric 
perturbations. However, it is believed that thin rings with large enough angular momentum 
will suffer from a GL-like instability [TS]. Since this would involve breaking the rotational 
symmetry of the ring we cannot investigate it using our methods. 

Finally, we consider the "doubly spinning" black rings of Ref. [21]. These also can be 
classified as "thin" and "fat". Presumably the arguments of Refs. [191 120] could also be 
applied to these but this has not been done so our work is the first study of the stability of 
doubly spinning rings. As in the singly spinning case, we find that fat black rings are unstable 
and we find no evidence of a rotationally-symmetric instability of thin rings. 

Notation 

Initial data for the vacuum Einstein equation in d dimensions will be denoted (S, hab, Kab) 
where S is a c? — 1 manifold with Riemannian metric hab and extrinsic curvature Kah- We 
shall denote the apparent horizon on the initial data surface by S. Va denotes the connection 
associated to hah- When discussing stationary black hole or black string solutions, S will 
denote a surface of constant t (in coordinates adapted to the timelike Killing field), which 
passes through the bifurcation surface where the past and future event horizons intersect. 
We will often consider 1-parameter families of initial data. The parameter will be denoted A 
and a derivative with respect to A denoted by a dot. For example, hab denotes a linearized 
perturbation of the metric of the initial data. An overbar denotes a quantity defined with 
respect to the unperturbed solution. We will work in units of G = c = 1. 

2 Black string instability 

Consider the black string in(i = n + 3>4 dimensions in the standard Schwarzschild 
coordinates: 

ds' = -f{r)dt' + ^rj-^ + r^dn'\^^ + dx^, /(r) = 1 - ^ , x~x + 27rL, (4) 



-i\r)de + ^^+r^dnl, 



where r = r_|_ denotes the location of the horizon. Consider a t = constant surface in this 
spacetime. To extend this surface to an Einstein-Rosen bridge with two asymptotically flat 
regions we define a new radial coordinate y such that 



1-1/2' 
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0<y<l. (5) 



In terms of this new coordinate the metric on our surface is 

We can now analytically continue y so that y G (—1, 1). In these new coordinates the bifur- 
cation surface is at y = and the two asymptotically fiat regions are y — )■ +1 and ?/ — )■ — 1 
respectively, y — > —y is an isometry which interchanges these regions. We shall denote this 
surface by S, with metric hab- The extrinsic curvature of S vanishes so (S, hah) provides 
time-symmetric initial data for the vacuum Einstein equation. 

We construct new time-symmetric initial data (S, hab) by the well-known method of con- 
formal rescaling. Let 

hab = ^^ hab (7) 

The Hamiltonian constraint reduces to Laplace's equation: 

V^^ = (8) 

where V^ is the Laplacian defined using hab-, and the momentum constraint is automatically 
satisfied. We seek solutions in which we excite just the lowest harmonic around the KK circle: 

^ = 1 + \f[y) cos{x/L) (9) 

where A is a parameter. The Laplace equation reduces to 

(l-yV(y)+4(n-%,(y) Arl ^ 

^^^'^ 2giy){l-y^) ^ ^^^ L^ g{y){l - y^)^ ^^^^ ^ ^ 

Solutions of this equation behave as e^''+/'''^^^~^ ^^ as y — )■ 1 and are smooth at y = 0. We 
fix f{%)) at some y = ymin < on the "other" side of the Einstein- Rosen bridge! so that 
fiUmm) = 1, and we choose f{y) to be the solution that decays as e~^'+^^^^^~y ^■' as y — > 1. We 
determine the solution numerically. Using this procedure we have constructed a 1-parameter 
family of initial data (E, hab) which reduces to the Schwarzschild black string initial data when 
A = 0. Note that {T,,hab) has the same asymptotic behaviour as (T,,hab) as ?/ — )■ 1. They 
differ in the other asymptotic region but this lies behind the horizon so this is not a problem. 
For small A, our initial data describes a small perturbation of the black string. 

Now we consider the local Penrose inequality ([3]). If the string is stable then our initial data 
must satisfy this inequality for sufficiently small A. Since our initial data is time-symmetric, 
the apparent horizon is an outermost minimal surface, i.e., it locally extremizes the area. 
Ref. |22] showed that the minimal surface is stable, i.e., the extremum is a local minimum|j 
Therefore it has area no greater than any surface that encloses it and hence Amin = v4app, the 
area of the apparent horizon. Furthermore, time-symmetry implies that the ADM energy is 
the same as the ADM mass, and so (IS]) reduces to 



^app < -^ (ii; 



^To obtain the data presented in this section we chose i/min = —0.1. We have checked that changing the 
location of this inner surface does not aher our results. 



^This point is discussed in more detail in Appendix Rl 



where M is the ADM mass of the initial data. Since f{y) decays exponentially, it follows that 
M is the same as the ADM mass of the unperturbed black string (i.e. M does not depend 
on A). The RHS of fITT]) is just the area of the horizon of this string and so we can rewrite 
this as /lapp(A) < Aapp(O) for sufficiently small A. Therefore we can prove that the string is 
unstable by showing that ylapp(A) > y4app(0) for arbitrarily small A. 

Since we are interested only in small A, we can expand Aapp = Aapp(O) + A/lapp(0) + 
(l/2)A^Aapp(0) + . . . where a dot denotes a derivative with respect to A. The first law of black 
hole mechanics applies to arbitrary linear perturbations that preserve the constraints on S 
|14] . The first law gives ^4(0) oc M(0). Here A can denote the area of either the event horizon, 
or apparent horizon, since they agree to linear order o For our initial data, M(0) =0 and 
hence 74app(0) = 0. Therefore our condition for instability reduces to 

iapp(O) > 0. (12) 

In Apppendix [B], we explain how to calculate ylapp(O). The result is: 



iapp(0)-^^^^^^Aapp(0) 



^^ ' {n + 2){n{n-l) + 2rl/L^) \dy , ^^^ 



(13) 



Hence we deduce that the black string is unstable if there exists a solution of the ODE (fTOl) 
that decays as ?/ — )■ 1 and satisfies 

, /I dfV n + 2 f , , 2rV\ , , 

We have solved f lTO|) numerically and Figure [1] shows ^ as a function of r_|_/L for the d = 5 
case. Note that A is positive for large r^/ L, consistent with stability, but becomes negative at 
small r^l L, indicating instability. In terms of r_|_, the critical value is r^/L = 0.8745. This can 
be compared with the critical value below which the GL instability exists|j Tj^/ L = 0.8762. 
The two results agree to an accuracy of 0.2%. 

We have performed an analogous calculation for Schwarzschild black strings with d = 
6, . . . , 11. The results are shown in Tabled! Our value is always smaller than the GL value 
and the agreement between the two values becomes better as d increases. 

These results are rather surprising. Violation of the local Penrose inequality is a sufficient 
condition for instability, but not a necessary one. For any given choice of initial data, the value 
of r^/L at which A changes sign should be smaller than the value at which the GL instability 
appears, as we find. Only by considering a sufficiently general class of initial data would one 
expect our approach to be able to identify the critical value of r^/L exactly. Nevertheless, the 
very simple initial data we constructed using the conformal rescaling method gives a value of 



^ To see this, note that two effects might contribute to the change in A at hnear order. One is the change 
in the metric at y = 0. The other is the change in the position of the horizon. But since y = is a minimal 
surface, the latter effect appears only at second order. Hence one has just the first effect, which does not 
depend on which kind of horizon we are discussing. 

^ We could not find a result to this accuracy in the literature so we determined it ourselves. The same 
applies for the other values quoted in Table [TJ 
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Figure 1: ^ vs. r+/L for the d = 5 black string. For other dimensions the plots look 
quahtatively similar. For large values of r^/L, A is positive but it becomes negative at a 
certain critical value, signalling an instability. For any number of dimensions this critical 
value is always smaller from critical value for the onset of the GL instability. In d = 5 we 
find (r_|_/L)crit = 0.8745, which differs by less than 0.2% from the GL critical value. The 
agreement gets better as the number of dimensions increases. 



r_|_/L which is within 0.2% of the critical value. This suggests that our initial data is (at first 
order) quite close to initial data for the GL unstable mode with the same x-dependence as our 
perturbation^ However, the GL unstable mode vanishes at the bifurcation surface whereas 
our initial data is non-vanishing there. Nevertheless, both perturbations are localized near 
the horizon so perhaps this is the reason for the surprising accuracy of our result. 

Note that we needed to work to second order in perturbation theory to demonstrate 
instability even though the GL instability is present at first order. This will be true also in 
the other examples we study. The reason is that our argument involves the change in mass 
and horizon area sourced by the perturbation. In general, a first order perturbation can lead 
to a first order change in the mass, angular momentum, horizon area etc. but these changes 
are governed by the first law [13]. This implies that the local Penrose inequality is saturated at 



first order. Given a first order metric perturbation bgah-, let ^g^-^ be a "non-dynamical", i.e., 
time- independent, perturbation that is obtained by a first order variation of the parameters 
of the unperturbed black hole solution, and has the same mass and angular momenta as bgah 
at first order. Write bgah = ^Qab^ + ^9ab where 5g^^ is "dynamical", i.e., time- dependent. If 
an instability is present then it is associated to 5g^^, not Sg^^,^. But, by construction, 5(7^ 



^More precisely, one can superpose the GL unstable mode, proportional to e^*, with its image under the 
i — > — t isonictry to obtain a time-symmetric unstable mode, which is more natural to compare with our initial 
data. 



d 


5 


6 


7 


8 


9 


10 


11 


(r+/L)crit 


0.8745 


1.2665 


1.5779 


1.8454 


2.0837 


2.3006 


2.5007 


{r+/L)GL 


0.8762 


1.2689 


1.5808 


1.8486 


2.0872 


2.3041 


2.5044 



Table 1: Critical value of r+/L obtained with our method (second row) compared to critical 
value that signals the onset of the GL instability (third row) as a function of the total number 
of spacetime dimensions d. The agreement improves as d increases. 

makes only a second order contribution to the mass and angular momentum. Hence we have 
to work to second order in order to see the instabihty using a local Penrose inequahty, even 
though this instabihty is present in the linearized theory. 



3 Rotating black holes 

3.1 Penrose inequality 

Myers-Perry black holes are uniquely parameterized by their mass M and angular momenta 
Ji where i = 1, . . . ,N = [[d — l)/2]. Black rings can be parameterized uniquely by M and 
a pair of dimensionless parameters (i^, a), which we shall define below. These quantities 
determine the angular momenta Ji, J2. However, (M, Ji, J2) do not uniquely specify the ring. 
If we define 



A = det 



-9(M,Ji,J2) 



diJuJ2 



then we can divide rings into "thin" rings with A > and "fat" rings with A < 0, each of 
which is uniquely parameterized by M and Jj. In discussing stability we can regard these as 
distinct families of solutions, e.g., a stable thin ring must remain on the thin ring branch if 
perturbed. Let Abi{{M, Jj) denote the area of the event horizon for the particular family of 
black holes under consideration. 

Consider a black hole belonging to one of these families and initial data corresponding to 
a small perturbation of the black hole. If the black hole is stable then, under time evolution, 
the perturbation should disperse, and the spacetime will settle down to a black hole belonging 
to the same family, with a small change in parameters. Let Mp and Jip denote the final mass 
and angular momenta, and Aj denote the area of the intersection of the event horizon with 
the initial data surface. As before, we define Amm to be the greatest lower bound on the area 
of any surface that encloses the apparent horizon on the initial data surface. Then 



^min <Ai< Abr{Mf, Jrp) < ^Bh(^, Jrp) 



(16) 



where E is the ADM energy of the initial data. The first inequality applies because the event 
horizon lies outside the apparent horizon, the second inequality is the second law and the 
final inequality uses Mp < E and the fact that Abh is an increasing function of mass at fixed 
angular momenta (which follows from the first law). 



The inequality ( !T6|) is not very useful because we know nothing about Jip. However, Ref. 
observed that one can circumvent this problem by imposing symmetries on the initial 
data which ensure that angular momentum is conserved. 

In our case, the families of solutions under consideration admit N commuting rotational 
symmetries, for which the angular momenta Ji are given by the associated Komar integrals. 
If we assume that our initial data preserves these rotational symmetries that these angular 
momenta will be conserved. Hence, for initial data corresponding to a small perturbation of 
the black hole, preserving the A^ rotational symmetries, stability implies the local Penrose 
inequality 

Amin < Abh(^, J^) (17) 

where Ji are the angular momenta of the initial data. 

This is saturated by a constant t slice through the unperturbed black hole. Hence a 
stable black hole is a local maximum of horizon area at fixed mass and angular momentum, 
and a local minimum of mass at fixed horizon area and angular momentum, in the space of 
rotationally symmetric, asymptotically flat, initial data. 

Finally, we have the problem that A^^^ is difficult to calculate. This is a problem that can 
be overcome by imposing an additional discrete symmetry on the initial data that ensures that 
the apparent horizon is a minimal surface. In 4d, an appropriate symmetry is "t— symmetry" 
[T^ 123] • This means: the initial data (Jl^hahiKab) is axisymmetric and one can introduce 
coordinates on the initial surface such that (i) the Killing field associated to axisymmetry is 
$° = {d/d(j))°', and (ii) — ?■ —0 is a diffeomorphism which preserves hab but reverses the sign 
of Kab- A surface of constant t in the Kerr spacetime has this symmetry. The existence of 
this symmetry implies that 

Kab = 2J(a$6) (18) 

where Ja is axisymmetric with Ja^"" = (and hence K = 0, i.e., the slice is maximal). 

The analogous symmetry for d > A dimensions, which we shall call "t — (pi symmetry", 
assumes TV rotational symmetries, generated by <!>*" = {d/d(j)i)°', and that the diffeomorphism 
(pi — > —(pi (simultaneously for all i) preserve hab but reverses the sign of Kab- Surfaces of 
constant t in all of our black hole families possess this symmetry. This symmetry implies that 

Kab = 2Jla%) (19) 

where J* is invariant under the rotational symmetries, and orthogonal to $^ for all j. 

Let S denote the apparent horizon in initial data with this symmetry and let n"" denote 
the outward unit normal to S in E. The condition that S is marginally outer trapped is 

K + {h^'' - n^n'')Kab = 0, (20) 

where K is the trace of the extrinsic curvature of S regarded as a surface in S. Now S must 
be invariant under the rotational symmetries hence n°'^a = 0- Using ( !T9|) we then deduce 
that the apparent horizon satisfies K = 0, i.e., it is a minimal surface. 

Next we must show that the apparent horizon is a stable minimal surface, i.e., a local 
minimum of the area rather than just an extremum. For d > A dimensions this has been 
established only for time-symmetric initial data [22]. However, for t — (pi symmetric initial 
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data, we prove in Appendix |X] that S* is a local minimum in the set of all homologous surfaces 
which lie outside S and are tangent to the Killing fields K^. Since the event horizon is such 
a surface, we can deduce that Aapp (the apparent horizon area) is a lower bound for Aj and 
hence 

Aapp<ABH(^,Ji). (21) 

To summarize: if the black hole is stable then this inequality must be satisfied by asymptoti- 
cally flat initial data describing a small perturbation of the black hole which preserves t — (pi 
symmetry. E is the ADM energy, and Jj are the angular momenta, of the initial data. The 
function A^n is defined by the family of stationary black holes that one is considering. 

3.2 Initial data 

We will construct initial data describing a t — 0j symmetric perturbation of the black hole 
by the Lichnerowicz method of conformal rescaling. We take the t — (j)i symmetric initial data 
(S, hab, Kab) on a constant t slice of our black hole solution and rescale 

hab = ■^''^"-'■^hab, Kab = "^-^K^b (22) 

where \E' is independent of the angles 0j. The momentum constraint is automatically satisfied. 
The Hamiltonian constraint reduces to 

V^^ - ^4^^R (* - ^-3-4/(d-3)^ ^ Q (23) 

where V is the connection, and R the Ricci scalar, associated to hab- 

Let U denote a surface homologous to the apparent horizon of the unperturbed black hole 
initial data, and lying a finite distance behind the horizon. Let S' denote the region of S 
exterior to f/, so 9S' = f/ U 6*00 where S^o denotes a sphere at spatial infinity. We will solve 
the above equation on S'. For small enough initial data, the apparent horizon will be close 
to that of the unperturbed black hole and therefore it will lie outside U and hence on T! . 

Now we consider boundary conditions. We will demand that \E' = 1 + 0{r~^^'^) at 6*00, 
which ensures asymptotic flatness. On U we will impose Dirichlet boundary conditions, i.e., 
we will specify '^\u- We will choose this function according to the kind of instability that 
we expect. The Hamiltonian constraint for the unperturbed black hole implies that -R > 0, 
which guarantees that there will be at most one solution of the above equation satisfying 
these boundary conditions. 

It would be straightforward to solve this nonlinear problem numerically. However, since 
we are interested in arbitrarily small initial data, we will proceed perturbatively. We seek a 
1-parameter family of solutions ^{X) with ^1^(0) = 1 (we suppress the dependence of ^ on 
the coordinates). We then Taylor expand \i/(A) = 1 + A\i/(0) + (1/2)A^^(0) + . . . where a dot 
denotes a derivative with respect to A. Substituting into the above equation gives 

V^^ - M = 0, (24) 

and 

v2§ -m = - ('^A^\ ^ ^2 (25) 

\ a — 3 / 
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with evaluation at A = understood. We solve fl2^ subject to the boundary conditions 
4/ = 0(r"'^+^) at 5*00 and \i/ = 4/\u on U for suitably chosen ^|r/. We then solve ( 1251) with 
the boundary conditions ^ = (9(r~'^"'"'^) at 6*00 and ^ = on Uu Our problem therefore is 
specified by the choice of the function \l/|[/. 

The boundary conditions at 6*00 imply that our initial data has the same angular momen- 
tum J as the unperturbed black hole solution, i.e., J{X) = J- The ADM momentum vanishes, 
so E = M, the ADM mass. The mass M(A) is determined by M(0) = M (the mass of the 
unperturbed black hole) and 

M(0) = ^^^1^ (r^-^v^) 1.^00, M(0) = ^-^^^^ (r^-^^) |™ (26) 

The first law (proved for t—(j) symmetric data in [15] and general initial data in [2]) guarantees 
that the first order change in the apparent horizon area satisfies 



Ti,pp(0) = M(0) (27) 



where T is the temperature of the unperturbed black hole. The calculation of the second 
order change is described in Appendix |Bl Next we expand the RHS of ( 12T]) to second order 
in A using 

^,AMM{X),J)] = ^(M(0),J)M(0) + ^^(M(0),J)M(0)2 

= ^Af (0) - ^M(0)2 (28) 

where 

fdM\ , , 

is the heat capacity at constant angular momentum of the unperturbed black hole. The first 
order terms in the local Penrose inequality (J2T|) cancel using the first law. At second order it 
becomes (for small enough A) 

Q > (30) 

where 

Q ^ M(0) - ^Ti,pp(0) - ^m? (31) 

The Penrose inequality is violated for arbitrarily small A if Q < 0. Our strategy for demon- 
strating instabilities will be to seek a function ^\u such that Q < 0. 



^This could be generalized by taking \l/ to be some specified function on U . The resulting solution would 
differ from our solution by a function ^2 satisfying the same homogeneous equation as ^. However the 
contribution from ^2 to the Penrose inequality expanded to ©(A^) is equivalent to replacing ^ with ^ + 
(1/2)A*2 in the 0[X) terms. Since the latter quantity satisfies the same equation as "f, the first law will 
ensure that these terms drop out of the Penrose inequality (see main text). Hence including ^2 has no effect 
and there is no loss of generality in our choice of boundary condition for 'f. 
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Note the presence of cj in the denominator of the final term of fl3T|) . What happens if cj 
vanishes for some particular black hole solution? This does not happen for the Myers-Perry 
solution. However, some black rings do have vanishing cj. In fact, "thin" rings have cj < 
and "fat" rings have Cj > 0, with cj passing through zero as one moves from the thin to the 
fat branch!^ For thin rings, the final term in ( 13T|) tends to make Q more positive. For fat 
rings it makes Q more negative. In particular, for a fat ring with very small cj it seems that 
Q will be negative as long as the initial data has M(0) 7^ 0. This strongly suggests that such 
a ring will be unstable. 

This is very similar to the argument used in Ref. [19] to deduce that fat black rings are 
unstable near the "turning point" where cj = 0. However, in place of the unspecified manifold 
Ai of Ref. |T9], we are working in the more well-defined setting of a 1-parameter family of 
initial data for Einstein's equation. This makes it clear that the instability (if one exists) is 
classical. Furthermore, the maximization of entropy argument of Ref. [12] is given a more 
precise formulation by the local Penrose inequality for rotationally symmetric initial datajj 

Finally, note that this argument is not a proof that fat black rings with small Cj are 
unstable because we have not excluded the possibility that other terms in (13T]) also become 
large. To demonstrate instability we will have to construct initial data for which Q < 0. 

4 Myers-Perry black holes 

4.1 Background solution 

We will consider a singly spinning Myers-Perry black hole with metric 



2 



ds^ = -^dt' + 5!^]^ (rf0 _ ndtf + ^dr' + p'de'' + r' cos^ ^rffij,_4) , (32) 



where 



A = r^ + a^ ^ , p"^ = r'^ + a^ cos^ 6 



j,d—5 



d-3 



2 /„2 I „2n2 „2a„- 2/1 n M 



E' = {r' + ay-a'Asm'9, Q 



r.,r a 



(33) 



^2^d-5 



The solution depends on the two parameters tm and a. The horizon of this black holes is 
located at r = r_|_: the largest root of A(r) = 0. For d = 4 and d = 5, the rotation parameter 
a have the bound a^ < r\,j/4 and a^ < r|^, respectively, with a strict inequality in the d = 5 
case since we want to consider only regular spacetimes. For d > Q, there is no bound for the 



^ A plot of temperature against mass at fixed angular momenta has a vertical asymptote at a point where 
cj = 0. Near the asymptote there will be two solutions with the same mass and angular momenta, these are 
the thin and fat rings. 

^Ref. [19] was uncertain about whether the "turning point" method requires rotational symmetry. We 
believe it to be essential for without it one cannot obtain a local Penrose inequality, which is the reason 
for believing that a stable black hole must be a local maximum of horizon area for fixed mass and angular 
momentum. 
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rotation parameter. The symmetry of this spacetime is Rt x U{1) x 5*0(^ — 3), where Rt is the 
time translation symmetry, U{1) is the rotational symmetry generated by d^ and SO{d — 3) 
is the symmetry of (ifi?^^-, part of the metric. The thermodynamical parameters are 

loTT a — 2 ^ ^ , . 

^ id- 3)rl + {d- 5)a' ^ ^^^ ^'^> 

47rr+(r^ + a^) ' r\ + a^ ' 

where M, J, A^, T and VLh are ADM mass. ADM angular momentum, area of the event 
horizon. Hawking temperature and angular velocity of horizon, respectively. ^1^-2 is the area 
of a unit [d — 2)-sphere. The heat capacity at constant angular momentum is 

l {d- 2)nd.2ri-\rl + a^f{{d - 3)r^ + [d - h)a^} 
^-^ 4 (rf-3)r| -6r2^a2 + 3(d-5)a4 ' ^ ' 

This quantity is negative for small a/rj^ and it never vanishes. However, for c? = 5 it diverges 
(and changes sign) at a/r+ = l/\/3. For o? = 6 it diverges at a/rj^ = 1 but it does not change 
sign. For ci > 7 it is finite and negative for all a/rj^. Note that a divergence (and sign change) 
also occurs for the Kerr black hole. 

Consider a surface of constant t in this geometry. To extend it through the Einstein- Rosen 
bridge, we define a new "radial" coordinate z by 

z' = '-—^ . (36) 



The resulting surface S has geometry 



ds' = '^ ^ d<p'' + ^^^^^dz"" + pHe^ + r^ cos^ ednl^_^) . (37) 



2_ ^^1 u^^2 , ^'+^ r ^^2 , „2j/32 , „2„„„2/3jr>2 

v2 ^ A/^^ ^^2 



Note that near the horizon we have A/z ~ A'(r+)r^ and therefore, above metric is regular 
at 2; = (the horizon). In these new coordinates the region of z < corresponds to the 
opposite side of the Einstein- Rosen bridge and the two asymptotically flat regions are related 
by the isometry z ^ —z. 

The extrinsic curvature of this surface has the form (fT9|) with just a single non- vanishing 
J*, corresponding to the rotational Killing vector $" = {d/d(j))°': 

Kab = 2J(,<l>b), Ja = y=dan (38) 

2pv A 
One can check that Ja is smooth at z = 0. 

4.2 Numerical strategy 

Now, we explain the numerical strategy to find the ultra-spinning instability of Myers- 
Perry black holes. First of all, we have to solve equations for the conformal factor Eqs.(l2 
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and fl25l). It is convenient to define new variables as 



From the asymptotic flatness, variables "$ and "^ must decay as 0{r~'^~^^) at inflnity and, thus, 
we impose Neumann boundary conditions for tp and ip at inflnity. Then, from asymptotic 
values of ip and ip, we can easily determine the deviation of mass using Eq. fl26|) . 

We also have to impose boundary conditions at an inner boundary. Without loss of 
generality we will choose the surface U deflned above to be located at 2; = —0.5, which is 
well behind the horizon {z = 0) of the unperturbed spacetime. Any smooth scalar in this 
geometry must be an even function of 6 near ^ = and an even function of 7r/2 — 6 near 
6 = tt/2. Such a function can be extended to an even function of 6 with period vr and can 
therefore be expanded in a Fourier cosine series with terms cos 2n6. Therefore we choose our 
boundary condition on U to be 

N 

^\u = l + J2^r.cos2ne , (40) 

where c„ (n = 1, 2, ■ ■ ■ , A^) are constants and the integer A^ represent a truncation of the 
Fourier expansion. We have flxed the normalization of ^, by choosing the flrst term in the 
expansion (l40l) to be unity. 

As explained above, we will set t(j\u = 0. We then solve equations (121|) - (125]) ) numeri- 
cally. The smoothness conditions just mentioned require that we impose Neumann boundary 
conditions at ^ = 0, 7r/2: 

deip\0=o = <9eV^|0=o = d0tp\0=^/2 = dBi'\e=T,/2 = . (41) 

Our task is to choose the constants A^ and c„ appropriately to flnd an instability. We shall 
take A^ = 4 and choose the constants c„, to minimize (numerically) the quantity Q deflned in 

Eq.daD. 

In order to suppress the truncation error at inflnity in our numerical calculation we use a 
compact radial coordinate w deflned as 

z = tan (-w) . (42) 

The coordinate range of w is tan^^(— 0.5) < w < 1. Discretizing Eqs. (l24|) and (125|) in the 
w-coordinate, we obtain systems of linear equations. We solved these equations using two 
different methods: conjugate gradient and successive over-relaxation [22]. We found that 
these methods gave same result within numerical error. 

In Figure. [21 as an example, we show the solution for d = Q, a/rM = 2Jlj and (ci, C2, C3, C4) = 
(3.51,3.12, 1.50,0.239) (These values of c„'s were chosen to minimize Q.). Since we have de- 



^'^ The background solution can be parameterized with rM and a/vM- The former is dimensionful and so it 
just defines a scale. Therefore we shall just state the value of the dimensionless parameter a/rM when giving 
results. 
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Figure 2: The solutions ip (left) and ip (right) for a/rM = 2 and (01,02,03,04) = 
(3.51, 3.12, 1.50, 0.239). The dotted line shows the position of the horizon of the unperturbed 
solution {w = 0). 



fined ip and 'ip as in Eq. (l39|) . these functions approach constant values at infinity w — )■ 1. From 
the asymptotic values of ip and ip, we can determine the deviation of mass M(0) and M(0) 
using Eq. (126 p . From Eq. (I94|) of Appendix [Bl we can obtain the first order deviation of the 
area of the apparent horizonO The explicit expression is 



^app(O) 



where we have defined 



and 



(ixJo(x)^(z 



X = cos 26 



0,x) 



{d - 2)7r(],_4rf "(r^ + a 
d-3 



1 +x 



(43) 

(44) 
(45) 



The first order changes in the mass and horizon area must satisfy the first law fl271) (recall 
that our perturbation does not change the angular momentum) and therefore we can use it as 
a global measure of the numerical error of our calculations. In the following we only present 
data for which the numerical error is less than 1% (see Appendix [C] for the details and the 
convergence tests). 

To evaluate the second order deviation of the area of the apparent horizon, we need to 
find the apparent horizon in the perturbed initial data. Let the first order position of the 
apparent horizon be 2; = XZ{x). From Eq. (1981) of Appendix [Bl we have the equation 



(46) 



d^Z dZ 

—— + A{x)^ + B{x)Z + C{x)d,i!{z = 0, x) = 
dx'^ dx 



Note that z plays the role of the coordinate r used in Appendix [ 
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where functions A, B and C are defined as 
{d-h)-{d-l)x 



A 
B 

C 



2i 2 

' \a X 



2(1 -x2) 

- 00 2/ 2 T 2N2n ^ K^ ~ 2 '^+ + ^ ~ ^^ J^ 

32r^(rj, + a^)^(l — a;"^) 

+ {4(rf - 2)rX + 7(rf - 3)r^a2 + (3rf - ll)a^}] 

2(rf-2)(d-3)r^ + (t/-5)a2 



(47) 



d-3 16ri(l-x2 



Regularity at x = ±1 gives the boundary condition that Z should be finite at x = ±1. Using 
the solution for ip shown in Fig. [21 we find the location of the apparent horizon shown in Fig. 
[3l From the solution of Zix), we can calculate the second order deviation of the area of the 




-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 

X 

Figure 3: The first order deviation of the apparent horizon for the solution of Fig. O 



apparent horizon. From Eq. (p6|) of Appendix |Bl we obtain the explicit expression for Aapp as 

2 



A 



app 



dx 



h{x) [2d,i>Z + i> 



d-3 J ^^Q \dx 



h{x)Z^ 



(45 



where 



lAx) = 47rfirf_4 



4ri-\rl + a^){l-z) /I + 



{d - 3)rl + {d - 5)a'^ \ 2 

d-5 



X 



I^[x) = 47rfirf_4 



rt^ fl + x 



[{{d - 3)rl + {d - 5)a^}a^x 



(49) 



16{rl + a2) V 2 

+ {4{d - 2)rl + 7{d - 3)rla^ + {3d - ll)a^}] 



and Io{x) was defined in Eq.( H5]) . Substituting the values for M, M, Aapp and A^pp into 
Eq. ipT]) . we can evaluate Q. 
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4.3 Results 

Using above procedure, we can evaluate Q for each choice of A^ and {c„,}. To find initial 
data which satisfies Q < 0, we seek to determine the {c„} which minimize the dimensionless 
quantity 

Q 



Q 



D^M 



D' 



7r/2 



^^^(^ = 0, 



(50) 



The factor of D^ is included because Q will scale as the square of the amplitude of \&. We 
have fixed the normalization of \1/ by equation (140|) but the c„ might become large in the 
minimization process. Dividing by D^ reduces the chance of the minimization algorithm 
running off to large c„. 

To carry out the minimization numerically we have used the downhill simplex method of 
Nelder and Mead [21] (see also [23]) setting A^ = 4 in Eq. (1101)0 In Fig. 11 we plot Q^in 
against a/rM for d = Q. We can see that the Qmin is negative for a/rM > 1.933, which proves 
that singly spinning 6d MP black holes with a/vM > 1.933 are unstable. 

Ref. [6] found that a time-independent mode indicating the onset of instability appears 
at a/rM = 1-572 so the instability should be present for a/rM > 1.572. This emphasizes 
that our approach gives a sufficient condition for instability but not a necessary condition. 
If we considered a more general class of initial data, which would involve going beyond the 
conformal rescaling method used above, then we should be able to get closer to the bound of 
Ref. [6|. However, it is striking that our very simple approach yields results that are so close 
to the exact onset of instability. 



0.05 



Q 



min 



-0.05 



-0.1 



"-1 1 



1.5 1.6 1.7 1.8 1.9 2 2.1 2.2 2.3 2/ 



a/rM 

Figure 4: The minimum value of Q is plotted against a/rM- Omin is negative for a/rM > 1.933. 



We have repeated the above calculation for other values of d. For d = 5 we find that 
Q is always positive. Therefore our results are consistent with the stability oi d = 5 single 
spinning Myers-Perry black holes against rotationally symmetric perturbations, in agreement 
with Refs. [HE]. (However, Ref. [16] found that there is an instability which breaks rotational 



We have checked that the result does not change much for larger N. 
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symmetry.) In Table [2], we summarize the our results for 6 < d < 11. We give the critical 
value a/vM beyond which we predict instability and the corresponding results of Dias et al [S]. 
Our approach becomes less good at identifying the onset of instability as d increases. So 
presumably our (first order) initial data looks less like initial data for an unstable mode for 
larger d. 



d 


6 


7 


8 


9 


10 


11 


ajru 


1.933 


2.380 


2.635 


2.803 


2.934 


3.048 


a/rM(Dias et al.) 


1.572 


1.714 


1.770 


1.792 


1.795 


1.798 



Table 2: Second row: critical value of ajru at which the local Penrose inequality is violated. 
Third row: onset of the instability found in Ref . [6] . 



5 Singly spinning black rings 

5.1 Background solution 

The metric of a singly spinning black ring 



in the coordinates of Ref. 



IS 



ds' 



F{V) 
F{x) 



1 +v 
dt - CR^^d(t)2 



ny) 



R' 



{x - y) 



;F{X) 



G{y) 
F{y)' 



dx^ G{x) 



G{y) • G{x) • F{x)' 



where 



F(o = 1 + Ae, G(o = (1 - r)(i + ^0, c = \\{\ - v) 



l + A 



(51) 



(52) 



1-A 

The constant R has dimensions of length and sets a scale for the solution. The dimensionless 
parameters A and v must lie in the range 

< z/ < A < 1 . (53) 

Absence of conical singularities fixes A and the periodicity of 0i, 02^ 

2v 



A 



0< 



>\^V2 



<27r 



v/r^A 



(54) 



1 + z/^ ' -..,..- ^_^ 

The ranges of the coordinates x and y are — 1 < x < 1 and —Ijv < y < —1. The event horizon 
is located at y = —l/u = y^- The mass, angular momentum, horizon area, temperature, and 
angular velocity are 



M 



A 



H 



2 {l-v){l + v'^) 
8V27i^u^R^ 

(1-Z/)(1 + 1/2)3/2 ' 



J 



1 + v 
v/2 V(l-^)(1 + ^') 

V 



OViT 



4V27ri?i/ 



^H = ^ 




Z/)(l + Z/2) 



(55) 



2 1 + 1/ 



The heat capacity at constant angular momentum is given by 

_ 12v/27r> - l/2)z.^v/rT^ 

A "thin" ring has < u < 1/2 and cj < 0. A "fat" ring has 1/2 < i^ < 1 and cj > 0. A ring 
with i^ = 1/2 is called "minimally spinning" because it has the minimum J for given M. 

Consider the induced metric on a constant t surface S. There is a coordinate singularity 
at the bifurcation surface y = yh- This can be eliminated by defining 

z^ = y-yh . (57) 

In the z coordinate, the induced metric hab on S can be written as 






G{y) 2 4:Z^dz^ dx^ G{x 



2 ^/ X r „, - -r , xU'v'i 



+ ^7r^ + 



F{y) ^' G{y) G{x) F{x) 



G^R^{l + yf 
F{x)F{y) 



d<Pi, (58) 



where y = y{z). Near 2; = we have G{y)/z'^ ^ G'{yh) and so this metric can be smoothly 
extended to a new asymptotically flat region with negative z. This is related to the original 
region by the isometry z — )■ —z. The coordinate range of z is then — Zmax < z < z^^y^ where 



2;max = V^l ^Vh (59) 

5.2 Numerical strategy for singly spinning black rings 

In the ring coordinates x and y (or z\ asymptotic infinity corresponds to a single point 
X — ;■ y — >■ — 1, which makes it difficult to impose boundary conditions there. Therefore we 
introduce new coordinates (ri,r2) defined by 

'"I = ^ Vr^ ' ^2 = R- TTT^ ' (60) 

X — Y [z) X — Y{z) 

where 

^^'^^ 7T7:::z ■ ^=^Vt^ '"' 

The function Y{z) is chosen to satisfy Y{z) ~ y for 2; — )• 2;max and l^(-s) — > —00 for 2; — )• — ^max- 
In these coordinates, spatial infinity corresponds to ri — )■ 00 or r2 — ?■ 00 and the asymptotic 
form of the metric is 

ds"^ ~ drl + rld4>i + drl + rld(()2 , (62) 

where (0i,02) = (1 — z^)/a/1 — A(0i, 02 ) so that the periodicity of 0i and 02 is 2 vr. InFigure[5l 
we depict the relation between (ri, r2) and (x, z) coordinates. Note that these new coordinates 
treat the two asymptotic regions asymmetrically, with asymptotic infinity in the negative z 
region corresponding to the point ri = R, r2 = 0. This is not a problem because we will not 
be solving any equations near this point. 
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Figure 5: Adapted coordinates (ri,r2) for i? = 1 and u = 0.6. The thick curve zjzm^y^ = 
corresponds to the event horizon of the background solution. The region inside the thick 
curve corresponds to the opposite side of the Einstein-Rosen bridge. We impose a Dirichlet 
boundary condition at z/^max = —0.5 which is located behind the horizon z/zmax = 0. 



Now we consider the equations for the conformal factor (1241) and fl25p . Requiring asymp- 
totic flatness implies that "^ and \E' must decay as 0(l/(r^ -|- r|)) at infinity. Thus, we define 



lb = — ^ ^ 



■■ rj + r^ + R^ ■■ 

lb = — ^ ^ . 



(63) 



and impose Neumann boundary conditions for ip and ip at infinity. Note that, in the above 
equation, we factorized rf + rl + R^ instead of rf + rg to make the transformation regular at 
the origin ri = r2 = 0. 

We choose our inner boundary U to be the surface z/zj^^^ = —0.5, where we impose the 
following Dirichlet boundary conditions: 



^1 



u 



^1 



u 



. 



(64) 



As we shall see in the next subsection, this simple boundary condition is sufficient to capture 
the instability of fat rings. At the axes ri = or r2 = 0, we impose Neumann boundary 
conditions as 

dr^tp\r^=0 = <9riV^|ri=0 = dr2lp\r2=0 = dr2i'\r2=0 = . (65) 

Furthermore, in our numerical calculation, we introduce compact coordinates Wi and W2 
defined as /^ x 

r2 — -R = ctan ( — tf;2 ) • (66) 



n 



ctan ( —wi 



IT 



We prepare grids at even intervals in the Wi and W2 coordinates. The constant c is chosen 
so that A^in/A^'aii ^ 0.025, where iVau represents total number of grid points and A^in is the 
number of grid points inside the horizon z < 0. The coordinate ranges are < U7i < 1 and 
2tan~^(— i?/c)/7r < u;2 < 1- Using the Wi and W2 coordinate, we solve Eqs.(j2^ and (1251) . 



20 



As an example, we show solutions for v = 0.6 in Figure [61 Note that tp, ip and Wi,W2 are 
dimensionless. It follows that our results depend only on the dimensionless parameter p, not 
on the scale R (or, equivalently, the mass M). 

From the asymptotic values of ip and ip, we can determine the deviation of mass M(0) 
and M(0) using Eq.(l26D. 





Figure 6: The solutions of ip (left) and ip (right) for u = 0.6. The dotted line is the position 
of the horizon of the unperturbed solution {z = 0). 

In the background solution, the apparent horizon is at z = 0. At first order this will 
be displaced to some new position z = XZ{x). This is found by solving equation (l98l) of 
Appendix IHl which takes the same form as equation fH6|) . where the functions A, B and C 
are given by 

2x + v{l + x^) 



A 
B 

C 



;i-x2)(z/a: + l) 



4A(l + A)(l-a;2)(z/x 
3(1 ^ 



— ^ \2\vx^ + 2Az/{(l + 3A)z/ - (1 + X)}x 
- (1 - Xy - (1 - 4A - hX^)u - 3A(1 + A)] 



(67) 



V 



Au{l -a;2)(z/x + l) ' 



Figd shows the position of the apparent horizon for the solution ip of Fig. [61 

The quantities Aapp and Aapp are calculated as explained in Appendix [Bl The results are 
given by eqs (HHl) and (H9l) with the functions /q, /i and I2 now being given by 



/o(x) 
hix) 



l2'n^u^uX{l - A2) 



R' 



ux] 



167rV^(l-x^)i?^ 
(l + z/)(l-i/)2(l + i/x) 

Att^u^R^ 



A(l-A2) 



(68) 



A 



[2Az/x^ + 2Az/{(l + 3A)z/ - (1 + A)}a; 



(1 - z/)2(l + zyx)4 y t/A(l + A) 

(1 - A)i^2 _ (1 _ 4^ _ 5^2)^ _ 3^^^ ^ ^^j 
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Figure 7: The first order deviation of tlie position of the apparent horizon for the solution 
of Fig. H 



Substituting our resuhs for M, M, 



Aapp and Aapp into Eq.(l3l]) 



we can evaluate Q. 



5.3 Results 

In FigJHl we plot Q against v where the quantity Q is defined in essentially the same way 
as for Myers-Perry black holes: 



Q 



Q 



D^M' 



D' 



dx 



1 Vl 



-.^{z = 0,x) 



x^ 



(69) 



The most striking feature of this plot is the apparent divergence at v = 1/2. This is the 
property anticipated in the discussion following eq. f l3T|) : Q diverges because cj vanishes 
at z/ = 1/2. Our results confirm the expectation that fat black rings with i^ ^ 1/2 have 
negative Q and therefore are unstable, as predicted by Ref. [I9]. Furthermore, we see that 
Q is negative at least up to ly ~ 0.95, beyond which the numerical error exceeds the criteria 
explained in Appendix O Hence all fat rings for which we have reliable results are unstable, 
in agreement with the prediction of Ref. [20] . 

For thin rings. Fig. |8] shows that Q is positive so there is no indication of instability. 
This is for the particular choice of initial data defined by equation ( 16^ . It might be the case 
that more complicated initial data leads to Q < 0. Therefore we have repeated the above 
calculation for initial data defined by 



N 



M 



U 



i+E 



Cr, cosmi 



(70) 



n=l 



where cos^ = x, and choosing the constants c„ by trying to minimize Q, as we did for the 
Myers-Perry black holes. We took A^ = 4 and considered values of u such that 0.05 < i' < 0.45. 
However, in all cases we found Q > so there is no indication of any rotationally symmetric 
instability of thin rings. 
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Figure 8: Q against v. Q is negative for v > 1/2. This proves the existence of an instabihty 
of fat black rings. 



6 Doubly spinning black rings 

6.1 Background solution 

We can study the instabihty of doubly spinning black rings in the same way as for the 
singly spinning black rings. The metric of the (balanced) doubly spinning black ring is given 
in Ref. [21]. It is labelled by a constant k with dimensions of length, and two dimensionless 
constants {a, z/j^ which lie in the range 



<a < 1 , 2^/a< u <1 + 



a . 



(71) 



In the second inequality, the lower limit corresponds to an extreme horizon and the upper 
limit to a naked singularity. When a = the solution reduces to the singly spinning ring. 

The solution is written using coordinates (t,x, y, 0i, 02)- The coordinates (0i,02) have 
canonical periodicity 27r and the same interpretation as for singly spinning rings (e.g., the S^ 
of the ring horizon is tangent to d/d(j)2)- The event horizon is located at 



y 



2a 



Vh 



(72) 



The mass, angular momenta, temperature and angular velocities are given in Ref. [21]. We 

The (0, V') coordinates of this reference correspond to our 



13 



Here (a, v) correspond to {v, A) of Ref. 

1,02)- 
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find that the heat capacity at constant angular momenta is 

24A;V(1 - ayh){l + i^yh){l + aylf{-2u^ - (1 + a)i^ + 1 + 10a + a^)} 



cj 



a{l + yh){l -a) 
X [l - (1 + 13a)yh + (2 + 31a + 8a^)yl + (-2 - 40a - lOa^ + 4a^)?/^ (73) 

+ 2a(2 - 28a + 5a^)yl + 2a(l + a)(5 + 21a + 5a^)yl + 2a'^{5 - 28a + 2a'^)yl 

- 2a^(-2 + 5a + 20a^ + a^)yl + a%8 + 31a + 2a^)yl - a\l3 + a)yl + a^yl°] "^ 



This vanishes at z/ = i^o where 



^(-l-a + v^(^^T9)(9^^TTy) = z/o • (74) 



cj is negative for u < uq and positive for u > uq. The Jacobian fllSp also vanishes at u = UiVl 
and so thin rings are those with u < uq and fat rings those with u > uq. 

As in the singly spinning case, we shall consider a surface S given by a constant t slice 
of this spacetime, with induced metric hat- We extend through the bifurcation surface at 
y = Uh on T, using the coordinate transformation z^ = y — yh, The coordinate range of the 
new coordinate z is taken to be — ^Jmax < z < z^ax so that z < corresponds to a new 
asymptotically flat region and where z^ax = V~^ ~ Vh- 

6.2 Results for doubly spinning black rings 

The numerical calculations for doubly spinning black rings are similar to those for singly 
spinning black rings and therefore we will only describe them very briefly. The coordinate 
system we use is same as Eq. (I6OI) . but the parameter R for doubly spinning black rings is 
defined as E? = 2k'^{l + a — z/)/(l — a). For this choice of R, the metric becomes explicitly 
asymptotically flat as in Eq. fl62|) . Introducing new variables {ip, ip) = (r^ + rj + -R^)A;~^(\1/, \1/) 
and compact coordinates f l66|) . we impose the boundary condition flM|) on the surface U given 
by z/z^i^x = —0.5, as before. 

Then, the flrst order deviation of the apparent horizon can be determined from the so- 
lution of -ipL^ Once we have determined (ip,^), we can extract the flrst and second order 
perturbations of the physical parameters, M, M, Aapp, Aapp and Q. 

In Fig. |9] we plot the sign of Q (which is what signals the existence of an instability) as a 
function of parameters {u, a) that specify the black ring (as before, the parameter k just sets 
a scale for Q). For increased clarity, we have parametrized u in terms of a new parameter /3 
defined as 

z/ = (1 + a)/3 + 2v^(l - /3) (75) 

where 1 + a and 2^/a are maximum and minimum values of u for a fixed a so P lies in the 
range < /3 < 1. The points • and x correspond to Q < and Q > 0, respectively. In this 
plot we have included only data which satisfies the numerical accuracy criteria discussed in 



^^See footnote [5] for an cxplanatfon of why this must be the case. 

^^ For doubly spinning black rings, we omit the explicit equations corresponding to Eqs. (|67p and 
because they are too long and unilluminating. 
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Appendix O The dotted curve corresponds to A = (or cj = 0). Hence the left and right 
hand sides of the dotted curve correspond to thin and fat rings respectively. We can see that 
fat doubly spinning black rings are unstable. As in the singly spinning case, this is to be 
expected near to the dotted curve for the reasons explained below (1311) . But our plot also 
demonstrates that all fat rings (at least those for which we have reliable data) are unstable. 
There is no sign of any rotationally symmetric instability of thin black rings. 
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Figure 9: The signature of Q for various a and /3. The points • and x represent Q < and 
Q > 0, respectively. Points to the left of the dotted curve correspond to thin rings, points to 
the right correspond to fat rings. 



6.3 Global Penrose inequality and minimally spinning rings 

So far, we have been discussing local Penrose inequalities, i.e., inequalities restricted to 
initial data describing a small perturbation of a given black hole. However, one can also 
obtain a global inequality for black rings. Consider asymptotically fiat initial data that, when 
evolved, "settles down" to a black ring solution with mass Mp and parameters {up, ap). Then 
the usual arguments give 



-4min < -4ring(Mi;', Up, ap) < A^ng^E , Up, ap) 



(76) 

where E is the ADM energy of the initial data. The second inequality holds because A^\ag is 
an increasing function of mass at fixed v^a. Now the black ring with greatest horizon area 
for given mass is the one with u = 1/2, a 
{3/16)^/3/7^E^/^. Hence the inequality 



and hence A^i^g{E, up, ap) < A-^i^glE, 1/2, 0) 



^min ^ ttjtiZ/ , 



-J- ^0.183 

16 V TT 



(77) 



must apply to any initial data which settles down to a black ring spacetimeJ^^I 



^^ This is much more restrictive than the "standard" Penrose inequality, in which the RHS is given by the 
area of a 5d Schwarzschild black hole of mass E, which has the same form with a* = (16/3)y^87r/3 « 15.4. 
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Initial data describing a ring with v = 1/2 and a = saturates the above inequahty. Such 
a ring is singly spinning and has A = 0: it is neither thin nor fat but sits between these 
two branches. It is "minimally spinning" in the sense that it is the singly spinning ring with 
smallest angular momentum for a given mass. The stability of a minimally spinning ring is not 
covered by the analysis of the previous sections (since it is neither thin nor fat). Refs. [271120] 
argued that such a ring must be unstable by considering the effect of dropping a particle 
with zero angular momentum into the ring. This argument neglects the backreaction of the 
particle. A more rigorous argument for instability can be obtained from the local version of 
the above Penrose inequality. 

Consider a 1-paramter family of asymptotically flat, initial data, with parameter A, such 
that for A = it reduces to a constant t slice of a minimally spinning ring with mass M , 
angular momentum J, angular velocity Q and temperature T. If the ring is stable then this 
initial data must satisfy the above inequality for sufficiently small A, i.e., 

«(A) - ^^ < a. (78) 

for small enough A. This is saturated at A = 0: a(0) = a*. The ffist correction arises from 
d(0) 

• /n\ ^min(O) 3Amin(0) • 

4 ( ^f ^ 7/ n\ 3E(0) 4 (2 



TM^I^ 



>)--'°))-^K5--"^) 



aJ l,m_,m^ (79) 



M3/2T \ M J 

In the ffist line we have used -E'(O) = M and in the second line we have used the ffist law and 
the Smarr relation [18]. We have checked that there exist perturbations for which d(0) 7^ 0: 
the methods of the previous sections can be used to construct a perturbation with J(0) = 0, 
i?(0) 7^ 0. If we choose A to have the same sign as d(0) then a(A) = a* + Ad(0) + . . . will violate 
the Penrose inequality for arbitrarily small A. Hence minimally spinning rings are unstable. 
Note that we were able to reach this conclusion using ffist order perturbation theory 
whereas previously we have always had to work to second order. This is related to a break- 
down in our observation at the end of section [2] that one can decompose a linearized metric 
perturbation into non-dynamical and dynamical parts with the former arising from a varia- 
tion of parameters in the unperturbed black hole, and the latter having vanishing mass and 
angular momentum at ffist order. For a minimally spinning ring, it is not possible to do this 
because a variation of parameters always gives a perturbation with d = 0, i.e., one cannot 
vary E and J independently this way. The best one can do is to split a perturbation into a 
non-dynamical part and a dynamical part, such that the latter has, at ffist order, vanishing 
energy but, in general, non- vanishing angular momentum. The fact that the dynamical part 
(which includes any instability) has non- vanishing angular momentum at first order explains 
why the Penrose inequality is non-trivial at ffist order. 
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We have shown that fat rings are unstable against rotationally symmetric perturbations. 
Therefore fat rings will admit rotationally symmetric linearized perturbations with exponen- 
tial time- dependence. Such modes should be analytic in v (since the background is) and 
presumably correspond to quasinormal modes of thin rings (since these appear to be stable 
against rotationally symmetric perturbations). Hence at z/ = 1/2, by continuity we cannot 
have an instability which grows exponentially with time: it must be sub-exponential. This 
is supported by the results of Ref . [2U] , where the minimally spinning ring corresponds to a 
point of inflection of the effective potential for radial perturbations. This suggests that the 
linearized instability should grow linearly with time. 

7 Discussion 

We have described how violation of a local Penrose inequality can be used to demonstrate 
the existence of certain types of black hole instability. In all cases, we constructed initial 
data by the Lichnerowicz conformal method. The simplicity of this approach is its main 
advantage. We did not have to derive and solve a large set of equations governing gravitational 
perturbations, or worry about gauge issues. Instead one just solves the single linear PDE (12^ 
to determine the first order perturbation, and (12 5 1) to determine the second order perturbation. 

A disadvantage of our approach is that in general it cannot be used to identify precisely the 
onset of instability. In the black string case, our results were surprisingly accurate. However, 
in the Myers-Perry example, we found that the local Penrose inequality is violated at a value 
of the spin parameter somewhat higher than that at which the ultraspinning instability is 
known to appear. This could be improved by considering more general initial data. 

We have seen how a local Penrose inequality can be used to make the turning point 
argument for instability (cf Ref. [19]) more rigorous. In this case, the approach does predict 
precisely when an instability should appear. 

The approach we have used here could be used in many other situations. The stability of 
non-uniform black strings [SU], or localized Kaluza-Klein black holes [HT] could be studied this 
way. One could also study Myers- Perry black holes with multiple non-zero angular momenta. 

We have considered only vacuum black holes. The same kinds of argument could be 
applied to charged black holes. In order to derive a local Penrose inequality, one must have 
a situation in which charge cannot leave the spacetime (and so the charge of the final state 
is the same as that of the initial data). For an asymptotically (locally) flat spacetime, this 
will be the case if there are no charged fields, or all charged fields are massive. The Maxwell 
field (or, more generally, p-form fields) must satisfy the Gauss law constraint. If initial data 
is constructed by the Lichnerowicz conformal method used above then it is straightforward 
to solve the Gauss law constraint by suitable conformal rescaling of the electric field of the 
initial data. 

Our methods are potentially more powerful for asymptotically anti-de Sitter spacetimes. 
In this case, the usual boundary conditions imply that angular momentum (as well as energy 
and charge) is conserved. Therefore one would not need to assume rotational symmetry to 
obtain a local Penrose inequality for rotating black holes. 

Instabilities of charged anti-de Sitter black holes have been the subject of much interest 
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recently [32]. It would be interesting to consider the various kinds of instabilities that can 
afflict such black holes using the methods discussed here. 

Finally, it would be interesting to relate the approach described here to the approach 
of linearized stability analysis. Can one prove that violation of the local Penrose inequality 
implies the existence of a linearized perturbation that grows exponentially with time? 
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A Apparent horizon for t — 0^ symmetric initial data 

We will show that, for t — (pi symmetric initial data, the apparent horizon 5" is a local 
minimum of area in the set of homologous surfaces outside 5* which are tangent to the Killing 
fields ¥''. 

Proof. We have shown that the apparent horizon is a minimal surface, i.e., extremizes 
the area. So we need to consider the second variation of the area. Let be a positive scalar 
on S and define the vector field V^ on S hj V^ = (pn"' where n"' is the outward unit normal 
to S in S. The map p — )■ exp{tV) then defines a 1-parameter deformation of S specified 
by the function 0. Since S is minimal, the first derivative of the area vanishes under this 
deformation. The second derivative is 

where Q denotes the expansion of the outward null normal to the deformed surface. Ref. [22] 
showed that 

— j = -A0 + 2X^D^<\> + (g + D^X'^ - X,X^) <P = O(0) (81) 

where A denotes the Laplacian on S, Da is the metric connection on 5*, indices are raised and 
lowered on S, Q is a certain scalar on S depending on the (instrinsic and extrinsic) curvature 
of S and X" is the vector field 

X" = {h"'' - n''n^)Ki,ji''. (82) 

In the time-symmetric case, X"" = and O is self-adjoint with respect to the obvious inner 
product. In this case one can argue that S must be a stable minimal surface, i.e., ( ISUj) is 
non-negative [22]. If X"" ^ then O is not self-adjoint but its eigenvalues have real part 
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bounded below, the eigenvalue Ai with the smallest real part is called the principal eigenvalue 
and must be real and non- negative |2S1 1221 • For the initial data considered in section 121 X"' 
is a linear combination of the Killing vectors 

X" = n^Jl^''' (83) 

We shall consider only deformations of the surface for which the Killing vectors $^ remain 
tangent to the surface. This is equivalent to considering only functions invariant with 
respect to these Killing fields. Hence we have X°-Da4> = 0. Restricted to the space of such 
functions, O is self-adjoint. Therefore the eigenvalues of the restriction are real. Since they 
cannot be less than Ai, they must be non-negative. Expanding in (1801) in eigenf unctions we 
deduce that (PA/dt^ > 0. 

Of course we really would like to establish the strict inequality (PA/dt^ > 0. This could be 
done if we knew Ai 7^ 0. To prove this we could show Ai 7^ for the unperturbed spacetime, 
then it will also hold in the perturbed spacetime. Or if Ai = in the unperturbed spacetime 
we would expect the perturbed spacetime to satisfy Ai 7^ 0. Either way it is clear that S 
should be a local minimum of area. 

B Apparent horizon area 

In this Appendix, we will explain how to determine the second order change in the area 
of the apparent horizon. Although many of the formulae here could be written in a covariant 
form, we will derive them in a coordinate system adapted to the symmetries of the black holes 
of interest. 

Consider {d — l)-dimensional initial data which is cohomogeneity-2, i.e., it depends non- 
trivially on at most two coordinates. We shall denote these coordinates as (r, x) and assume 
a < X < b. We assume that the initial data has a metric of the form 

hab = n{X,r,xYhab (84) 

where hab is initial data for the black hole whose stability we are investigating. Furthermore 
we assume fi{0,r,x) = 1, so hab = hab when A = 0. In the cases we are interested in, r will 
be a "radial" coordinate. For the black string, x will be the coordinate around the KK circle 
and for the rotating black holes, x will be a direction cosine. Note that hab depends on A only 
through the conformal factor Q. 

Now consider a surface r = r{x) in this initial data. Its area is given by a functional of 
the form (using a prime to denote a x-derivative) 

/"' 

v4[A,r(x),r'(x)] = / dx F{X,x,r{x),r'{x)). (85) 

J a 

Note that ([84]) imphes that 

F(A, X, r{x),r'{x)) = fi(A, x, r(x))'^"^F(x, r(x), r'(x)) (86) 
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where F{x,r{x),r'{x)) = F{0,x,r{x),r'{x)). In the cases we are interested in, the apparent 
horizon is a minimal surface, i.e., it extremizes A. Let r = R{X,x) denote this minimal 
surface, which must satisfy the Euler-Lagrange equation 



E{X, X, R{X, x), i?'(A, x), R"{X, x)) = 0, 



(87) 



where 



E{X, X, r(x), r'(x), r"(x)) 



d dF OF 

dx Or' dr 



The area of the apparent horizon is 

Aapp(A) = ^[A,i?(A,x),i?'(A,x)]. 
Using the Euler-Lagrange equation, we have 



KpW 



dxF{X,x,R{X,x),R'{X,x)) 



V'^'^ / r(x)=R{X,x) 



(89) 



(90) 



where we are using a dot to denote d/dX (which acts just on the first argument of F). The 
second term is a surface term coming from the boundaries of the x-integration. For the black 
string, we impose periodic boundary conditions so the surface term vanishes. For the rotating 
black holes, we use dF/dr' = Q'^~'^dF/dr'. Functions r(x) describing a smooth surface will 
satisfy dF/dr' = aX x = a,b. Hence we can neglect the surface term. 
Taking another derivative gives 



^app(A) = dx 



d r)F r)F 
F(A,x,i?(A,x),i?'(A,x))-|^^-^ 



R{X,x) 



r{x)=R{\,x) 



(91) 



where R{x) = R{0,x). A surface term can be neglected for the same reason as before. To 
simplify this expression, note that ( 186|) implies 

F{X, X, r(x), r'(x)) = {d - 2) ^|^^'''S^|| f(A, x, r(x), r'(x)) (92) 



f2(A, x,r(x)) 



and 



F(A, X, r(x), r'(x)) = {d - 2) (n{X, x, r(x)) + {d - 3)h{X, x,r{x)f) F{x, r(x), r'(x)) (93) 
The first of these implies 



^app(O) = {d-2) dx n{0, X, i?(x))F(x, R{x), R'{x)) (94) 

J a 

To evaluate the second term in (l9Ti) . we use ( l92l) and the equation of motion ( IHTl) to obtain 

{d-2 

^ d ( Vl\ f ,dF 



d dF dF 

dx dr' dr 



r{x)=R{\,x) 



d_ (n\ d£_d_n 

dx \ ^ I dr' dr \ Q, 



{d-2) 



dr \Q 
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dr' 



r(x)=R{X,x) 

_d_ foX dF^ 
dx \ Vt I dr' 



(95) 

r{x)=R{X,x) 



Substituting into fl9T]) and evaluating at A = now gives 



i^pp(O) = {d-2) dxl ii + (d-3)f2 



on 

dr 



\=0,r{x)=R{x) 



F{x,R{x),R'{x)) 



A=0 



r'^ - F 
dr' 



dx 



A=0 



dF 
dr' 



R{0,x) 



(96) 



r{x)=R{x) 



where R{x) = R{0,x). Note that this expression involves the metric perturbation to second 
order in A (via Cl) but we only need to determine the first order perturbation -R(0, x) to the 
position of the apparent horizon. 

Next we need to determine R{0, x). To do this, we substitute ( 186|) into the minimal surface 
equation (187|1 . The result is 



(97) 



n^-^Eix,ri^)yix)yix)) + 1 {n^~^) ^ + 1 {^^ (r'Z - ^ 



dx 



dr' dr 



dr' 



r=R{\,x) 

where E{x,r{x),r'{x)) = E{0,x,r{x),r'{x),r"{x)). Evaluating at C(A) gives 
d 



■jtE{x, R{\, x), R\X, x), R"{X, x)) 
dX 



J A=0 



+ (rf-2) 



dn 

dx 



A=0 



dF 

dr' 



dQ 

dr 



A=0 



r'^ - F 
dr' 



(98) 



r=R{x) 



The first term gives an expression linear in R{0,x) and its first and second x-derivatives. 
The second term is a source depending on the first order metric perturbation Q. Hence this 
equation is a second order linear ODE that, with suitable boundary conditions, determines 
R{0,x). 

We now apply this to the d = n + 3 dimensional black string (j6]), for which 



F{x,r{x),r'{x)) 



UJr,.r\ 



4:r\{y'{x)Y" 



[I - y{x] 



a\n 



g{y{x)){l-y{x)^y 



(99) 



2 

where Un denotes the area of a unit round n-sphere, fi = \1/™ and \l/ is given by equation (|9]). 
The horizon is at y = in the unperturbed spacetime hence R{x) = 0. Equation (p4l) gives 



2itL 



^app(O) oc / dxcos{x/L) = 
Jo 

so the first order change in Aapp vanishes, as discussed above. Equation (I5S1) reduces to 

n^-1 fdf 



(100) 



i?"(0,a;)-^^^^^i?(0,x) 



2nr^ \dy 



cos{x/L) 



(101) 



j/=0 
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The general solution of this equation is a sum of a particular integral proportional to cos{x/L) 
and an arbitrary linear combination of e^^V"("~^)/(''+^^). Since no such linear combination 
obeys the required periodicity a; ~ x + 27iL, the only acceptable solution is 

Hence this is the first order perturbation to the position of the apparent horizon. We now 
substitute this into (15^ to obtain the second order change in the apparent horizon area (IT5]1 . 

C Numerical errors 

In this appendix we present some of the convergence tests that we have performed in order 
to check our numerics. As discussed in the text, the main quantity that we have monitored 
to check numerical errors is the First Law, which we know must be satisfied in the continuum 
limit and therefore it provides a global measure of the numerical error. 

The mass of the background spacetime M (0) sets the scale of problem, but the actual 
"size" of the perturbation can be captured by the quantity 

D^= [ 4l\l, (103) 

Jh 

which is evaluated on the horizon of the unperturbed spacetime. In order to measure the 
numerical error we want to consider a dimensionless quantity for which the size of the per- 
turbation has been scaled out. Therefore, the following quantity 



DM{0) 



M(0)_,^'^^PP(°) 



4 



(104) 



should be a good measure of the numerical error and this is what we have used. We reject 
data for which e > 0.01. Furthermore, e is an estimate of the error in Q because the equation 
for \l/ is the same as that for \E' but with a source. Since we want to know whether Q is 
positive or negative, we reject data for which \Q\ < lOOe (this is only an issue for black rings, 
e.g. \Q\ becomes small as z/ — )■ or 1 for singly spinning rings). 

We have written two codes, one based on second order finite differencing approximation 
and the other based on a pseudospectral collocation approximation. In Figure [10] we present 
the results of the convergence tests that we have performed for the pseudospectral code for 
Myers-Perry black holes in all dimensions J^l As this figure shows, the error e defined above 
decreases exponentially with the number of grid points A^, as expected for smooth functions. 
It is interesting to note that ioi d > 9 our code could not find acceptable solutions for low 
resolutions, i.e., 10 grid points, but for higher resolutions the behaviour of our code seems to 
be independent of the number of spacetime dimensions. To carry out the tests we have chosen 



^^For singly spinning rings and doubly spinning rings the results are qualitatively similar and we do not 
show them here. 
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Figure 10: Convergence tests for MP black holes in all dimensions. The error decreases 
exponentially as a function of the number of grid points independently of the number of 
spacetime dimensions d. 

the number of grid points to be the same in all coordinate directions. Finally we note that 
we have also observed second order convergence for the finite differencing code. Therefore, 
we conclude that our codes exhibit the expected convergence to the continuum according to 
each differentiation scheme and hence we believe that the numerics are under control. 
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